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Abstract. Two new expanding X-ray rings were detected by the Swift XRT instrument during early follow-up observations of 
GRB 061019 and GRB 070129, increasing to 5 the number of dust scattering X-ray halos observed around GRBs. Although 
these two halos were particularly faint, a sensitive analysis can be performed that optimizes the method originally developed 
by Tiengo & Mereghetti (2006) to analyze dust scattering rings observed with XMM-Newton for the Swift satellite. In the case 
of GRB 061019, a known giant molecular cloud is identified as the one responsible for the scattering process, and its distance 
is accurately measured (d=940±40 pc) through the dynamics of the expanding ring. In the second case, XRT observed both 
the main peak of the prompt emission of GRB 070129 and the scattering halo, but the small number of detected halo photons 
prevents us from distinguish between different dust models. 
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1. Introduction 

The scattering of X-rays from interstellar dust grains is of- 
ten responsible for the formation of diffuse halos surround- 
ing bright sources. The study of such halos can then be used 
to infer the properties of the dust and its spatial distribution 
( Mat his & Lee 1991; Draine 2003 ). and, in case of variable 
sources, can also yield inform ation on the source's distance 
dTriimper & Schonfelderlll973h . 

Thanks to rapid follow-up observations with sensi- 
tive X-ray telescopes, it has recently become possible to 
observe halos around gamma-ray bursts. Scattering halos 
due to dust in our Galaxy have b een detected aroun d 
three gamma-ra y bursts: GRB 03120 3 dVaughanet al.ll2004 . 
(Vaughan et al 



GRB 050724 



2006), and" GRB 0507 13A 



( Tiengo & Mereghetti 20061) . Due to the short duration of the 
bursts and the relatively small thickness of the dust layers, such 
halos appear as expanding rings. For the dust distances D^,,,, 
typically involved here (hundreds of parsecs) the rings have an- 
gular radii 9 of a few arcminutes after several thousand seconds 
from burst time: 



6[arcsec] 



827 



DduApc] 
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The brightness of the halo, Ihalo, depends on the inten- 
sity of the GRB emission, \qrb, and on the scattering optical 
depth, t: \ H alo = Igrb (l-e~ T ) - r l CRB . The same dust is 
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responsible both for the X-ray halos and for the interstellar 
reddening, so r is usually inferred from measurements of the 
optical reddening Ay- However, different relations have been 
reported in the literature between Ay and r. Draine has devel- 



oped a detailed physical model of the dust (lDraind l2003) and 
extracted from this model the following relation between op- 
tical e xtinction and scatterin g optical depth in the soft X-ray 
range (Dr aine & Bond 2004): 

T( ^ .16A y (JLf . 
However, a different relation was empirically found 



Predehl & Schmitt ( 1995), that is: 



(2) 
by 



t(E) * 0.056A 



v (i!v) 



(3) 



,lkeV 

If both Ihalo and Iqrb are known or can be es timated, as in 
the case of GRB 050724 dVaughan et al.l 120061) . some infor- 



mation on the amount of dust and its properties can be de- 
rived. Otherwise, by knowing or making assumptions about 
the amount of dust and its properties, the measurement of 
the radiation scattered in the halo can provide some infor- 
mation on the GRB prompt emission at a few keV. This 



has b e en done, for example, for G RB 031203 dVaughan et al 
2004 iTiengo & Mereghetti! 120061) . Here we report on two 



new dust-scattering rings recently obs erved with Swift ar ound 
the gamma-ra y bursts GRB 06101 9 dRomano et al.l l2006) and 
GRB 070129 (I Vianello et ai1l2007l) . Our analysis is based on a 
metho d particularly convenient in the case of faint, expanding 
halos dTiengo & Mereghettill2006l) . which is briefly described 
in the next section. 
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2. Analysis method 

The method developed by Tiengo & MereghettH ( 2006 ) for de- 
tecting and analyzing expanding dust halos can be applied to 
any imaging X-ray detector. We applied it to dat a obtained with 



GRB « miff 



the X-Ray Telescope (XRT) onboard Swift (Bu rrows et al 
l2000h . 

For each count detected by the XRT, we computed the ar- 
rival time t, in seconds from the burst time, and its squared 
angular distance from the GRB position, 8 2 . Based on these 
quantities, we built a dynamical image, that is a 2-D histogram 
with t on the x-axis and Q 2 on the y-axis. In such an image an 
expanding ring forms a straight line with an angular coefficient 
depending on the dust distance, as can be seen from EqQ] In 
this way we can easily detect a dust-scattering ring and estimate 
the distance of the dust layer from the line slope. In order to de- 
termine the halo intensity and dust layer distance and thickness, 
we removed the GRB afterglow and any other point sources 
from the event list, and we computed the value D - 827 t/8 2 
for each detected event. In the ideal case (thin dust layer and 
perfect instrument), every halo event will have the same value 
of D, corresponding to the actual distance of the dust layer. In 
real cases, instead, the instrument point spread function (PSF) 
and the intrinsic width of the scattering dust layer produce a 
distribution of D centered on the dust layer distance. If we plot 
the distribution of the quantity D (dust distance distribution 
hereinafter), we obtain a single peak centered on the real dust 
layer distance, formed by the photons belonging to the ring, 
superimposed on a background distribution. 

If the instrument PSF is described well by a King function 
(as in the case of both XMM-EPIC and Swift-XRT), the dust 
ring peak is well-modeled by a Lorentzian function. When the 
observation is continuous (as typically in the case of XMM- 
EPIC observations), the background distribution is a power law, 
so the dust distance distribution is described by a Lorentzian 
plus a power law. Unfortunately, the Swift observations are al- 
ways discontinuous, owing to the low earth orbit of the satel- 
lite and to the particular observation plan adopted by XRT 0. 
Since in this case the background distribution cannot be easily 
described analytically, we developed a Montecarlo simulation 
that reproduces the contribution to the dust distance distribu- 
tion produced by the background counts, as well as the dust 
scattering ring, taking the actual observation plan into account. 
This simulation is also useful for understaning some other as- 
pects of the phenomenon better, such as the effects due to the 
instrumental PSF and effective area or due to the dust layer 
intrinsic width. Our simulation can generate a GRB prompt X- 
ray emission based on input parameters describing its timing 
and spectral characteristics; alternatively, it can use as input 
an event file obtained during real observations of the X-ray 
prompt emission. Then it models the scattering process using 
the layer distance and thickness as input. Our simulation can 
use either of the relations (j2]i or ([3). Taking the characteristics 
of the instrument (PSF and effective area), the background and 
the given observation plan into account, the program generates 
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1 Typically the observing plan of a GRB afterglow is composed of 
several individual exposures separated by time gaps. 



Fig. 1. Image (top panel) and dynamical image (bottom panel) 
of GRB 061019 in the 0.2-4 keV energy range. We have ex- 
cluded from the images a circle with a radius of 30" centered 
at the GRB position to remove the afterglow contribution. In 
the upper panel the cross marks the afterglow position. 



an output event file in a standard format allowing its analysis 
as that of a real observation. 



3. GRB 061019 

GRB 061019 was discovered with the Swift Burst Alert 
Telescope (BAT) and localized at Galac tic coordinates 
1=181 .74° and b=4.26° dGuidorzi et al.l2006allbh . The XRT and 
UVOT began their observations of the burst location 2800 s af- 
ter the trigger ti me To _ 04:19:06 U T. A dust scattering ring 
was discovered (Romano et al . 2006) in the first three orbits of 
XRT Photon Counting data (from- r +2800sto~ r + 17000 
s). 

Our analysis of these data has confirmed the detection of 
the ring. The scattering halo is very visible in the sky image 
and in the dynamical image shown in FigQ] These images, as 
well as all the analysis discussed below, refer to the 0.2^1 keV 
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energy range. The line formed by the halo events in the dynam- 
ical image is rather wide due to the intrinsic thickness of the 
expanding ring, which remains clearly visible until ~ 23000 s 
after the burst. We estimated the background count rate from 
four regions of the XRT field of view, at least 4' away from the 
afterglow position and with no evident point sources. The aver- 
age count rate was (3.4 + 0.3) x 10~ 8 counts s _1 arcsec~ 2 and no 
significant spatial variability was found by comparing the dif- 
ferent regions. Using the Montecarlo simulator described above 
with this background rate (and its error), we estimated the dust 
distance distribution expected for the background counts. This 
distribution was subtracted from the data, taking both the sta- 
tistical (poissonian) errors and that of the background rate into 
account, in order to derive the net distance distribution for the 
halo shown in Fig. [2] After removing the background, the halo 
counts are ~ 180 ± 20 in the 0.2 - 4 keV energy range. The 
peak is well-fitted with a Lorentzian centered at D — 941^12 pc 
and with FWHM width AD = 427+^ 5 pc OrVd.o.f. = 0.74 with 
16 d.o.f.). Our simulations show that in this case the XRT PSF 
produces a peak width of ~ 180 pc. By fixing the Lorentzian 
FWHM at this value, we obtained a fit to the peak with the 
same center (D = 938 + 12 pc), but with a significantly worse 
chi-square (^ 2 /d.o.f. = 1.71 with 17 d.o.f.) corresponing to 
a chance probability of 4%, so that the hypothesis of a zero- 
thickness cloud is rejected at 3 sigma level. We thus conclude 
that the PSF is only partially responsible for the observed ring 
width and that the dust layer causing the scattering has an in- 
trinsic width ADtntr > 150 pc. 

The halo spectrum can be fitted with an absorbed power 
law. Due to the small number of counts we fixed the hydro- 
gen column density Nh to the best-fit value obtained for th e 
afterglow (N H = 0.9 x 10 22 cm 2 , iGuidorzi et ail l2006bl) . 
finding a photon index r = 1.8 ± 0.9 and a halo fiuence of 



~ 3 x 10 } cmr 1 Or /d.o.f. = 1.3 with 2 d.o.f). Adopting the 
optical depth-energy relation in Eq|2]and assuming the optical 
extinction reported in lSchlegel et alT (ll998l) for the bust direc- 
tion (Ay = 3.4), we estimated the prompt X-ray fiuence f ol- 
lowing the method described in lTiengo & Mereg hetti ( 2006J) as 
~ 2 x 1 0~ 7 erg cirT 2 (1-2 keV) . If we extrapolate the BAT spec- 
trum (r = 1.85+0.26, 15- 150 k eV fiuence = (2.4 ± 0.3) x 10~ 6 
erg cm JSakamoto et al . 2006) to the 1—2 keV energy range, 
we obtain a fiuence of (2.5+0.3) x 10~ 7 erg crrT 2 , in agreement 



with our estimate. 

We thus conclude that, for this GRB, the prompt emission at 
soft X-ray energies (0.3 - 4 keV) is consistent with an extrap- 
olation of the hard X-ray emission (15-150 keV), suggesting 
that it probably has the same origin. If we instead adopt the op- 
tical depth-energy relation in Eqf3]we obtain a flux in the 1-2 
keV band that is ~ 3 times greater than the BAT-extrapolated 
value. However, owing to the errors in the estimate of the spec- 
tral parameters of the ring and in the optical extinction Ay, this 
result is not conclusive. 

4. GRB 070129 

On 2007-01-29 the BAT instrument was triggered at 23:35:10 
UT by a precursor of GRB 070129, and the Swift satellite 
automatically slewed to the derived sky position. This al- 
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Fig. 2. Background-subtracted dust distance distribution of 
GRB 061019. The peak is due to the expanding ring. 

lowed detection of the main peak of emission, which was de- 
tected at 23:41:30 UT by both the BAT and XRT instruments. 
Despite the burst location at high Galactic latitude (1=157.2°, 
b=-44.7°), in a re gion of moderate ab sorption (optical ex- 
tinction Ay =: 0.4, ISchlegel et al. 1 19981) . a detailed analysis 
of the XRT data led to the discovery of a faint partial ring 
(IVianello et al.1 120071) centered at the afterglow position. The 
ring radius increased from ~ V to ~ 3' during the second and 
third Swift orbits (from ~ T + 4450 s to ~ T + 13100 s). 
Figure [3] shows the dynamical image obtained from the XRT 
Photon Counting data. The expanding ring is evident in the sec- 
ond and third orbits. The cross section of the scattering process 
decreases rapidly with energy (see Eqs|2]and[3]l and becomes 
vanishingly small above a few keV, so just as for the previous 
burst, we restricted our analysis to energies below 4 keV. From 
an annulus with radii of 5' and 8' centered on the afterglow and 
excluding point sources and flickering pixels, we estimated a 
background count rate of (8.9 + 0.7) x 10~ 8 cts s arcsec~ 2 . 
Inspecting the dynamical image, we selected the region con- 
taining the halo as an annulus of inner radius of 90" and outer 
radius of 160". If we consider for the moment only the second 
orbit, this region contains 45 + 6 counts, of which 11+3 are 
expected to be background. Thus the net number of halo counts 
in the second orbit is 32 + 8, corresponding to a detection at the 
~ 4<x confidence level. 

Owing to the small number of halo counts, we based our 
analysis on the integral distribution of D. This approach has the 
advantage of using all the available information and avoids the 
arbitrariness of data binning. The integral distributions of the 
Lorentzian peaks are arctan profiles, with Ddust corresponding 
to an inflection point. Consider the expected halo region, which 
is an annulus centered on the GRB position with inner radius 
of ~ 30" (in order to exclude the afterglow) and outer radius 
of ~ 300" (the scattering cross section is negligible at larger 
radii): the integral dust distance distribution obtained from such 
a region using only the second and the third orbits is reported 
in Fig. [4] In order to model this distribution, we used the XRT 
observation of the prompt X-ray emission as input for our sim- 
ulation, initially with the optical depth-energy relation in Eqf2] 
and the background rate estimated above. Then we defined a 
grid of possible values for the distance of the dust layer, and for 



4 



G.Vianello, A. Tiengo & S. Mereghetti: GRB halos 



Dynamical Image of GRB 070119 
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Fig. 3. Dynamical image of GRB 070129. The halo is visible 
in the second and third orbits. Is not visible in the first one due 
to the high background level and short exposure time. 

Integral distance distribution of GRB 070129 
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Fig. 4. Integral dust distance distribution of GRB 070129 (light 
gray points). The black line is the simulated distribution with 
the two scattering rings at D — 150 and D = 290. The dark gray 
line is the distribution expected from a uniform background. 

every value of the grid we made a large number of simulation 
runs in order to derive the expected shape for the integral dis- 
tribution. We then used the Kolmogorov-Smirnov (KS) test to 
compare the observed distribution with those generated by the 
simulation. We found that the model that described the obser- 
vations (with a KS-test probability of 55%) best corresponds to 
a dust distance of D — 290 pc. In contrast, if we try to describe 
the observed distribution with what is expected only from the 
background, we obtain a probability of only 5 x 10~ 5 (see 
Fig. 0}. We applied the same procedure to a late observation 
of the afterglow of GRB 070129, when the dust ring has disap- 
peared. In this case the distance distribution measured from the 
same annular region is described well (KS-test probability of 
99%) by the one expected from uniform background, demon- 
strating the accuracy of the background distributions derived 
with our simulation. Inspecting Fig. |H one can see that there 
is another inflection point at ~150 pc. This could be due to an- 
other, fainter dust ring. Using the same procedure as described 
above, we found a significant improvement (KS-test probabil- 
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Fig. 5. Image of the afterglow of GRB 070129 obtained with 
XRT during the second Swift orbit since trigger time. The im- 
age has been rebinned and smoothed to show the faint diffuse 
structures better. The dust scattering ring is visible around the 
afterglow. The ring appears asymmetric, since is brighter in the 
region north west of the marked line. 

ity of 90%) with two layers of dust, at distances D = 150 pc and 
D = 290 pc. We thus conclude that the halo observed around 
GRB 070129 is probably due to two dust layers along its line 
of sight. The knowledge of the prompt X-ray emission should 
allow us to test and possibly distinguish between different pro- 
posed optical depth-energy relations for the scattering process 
like Eqsf2]and Eqf3] corresponding to different dust properties, 
such as grain size distribution and composition. Our simula- 
tions with the best model derived above for the dust distribu- 
tions predicts 40 + 7 halo counts in the second XRT orbit us- 
ing the Draine optical depth-energy relation (Eq.|2} and 18+4 
halo counts using the Predehl relation (Eqf3j). Although the first 
value is closer to the observed one (32 + 8 cts), the inadequate 
statistics do not allow us to derive a firm conclusion in favor of 
one of the two possibilities. 

There is some evidence that the ring spatial distribution is 
not uniform around the GRB position, being more intense in 
the north west part. Dividing the halo region into two parts, 
as indicated by the line in Fig. [5] we find 24 + 5 counts 
(halo+background) in the north west part and 10 + 3 on the 
opposite side. We checked that this uneven distribution, corre- 
sponding to a chance probability of ~1.6%, is not due to a non 
uniformity in the XRT exposure. The asymmetry in the ring in- 
tensity is thus probably due to a inhomogeneous distribution of 
the dust in the distant layer. 

5. Discussion 



The f irst dust scattering ring around a GRB dVaughan et al 



2004) was disco vered thanks to the rap id localization provided 
by INTEGRAL dMereghetti et al.ll2003l) and the large collect- 
ing area at the X-ray energy of the EPIC instrument on XMM- 
Newton. The results reported here demonstrate that, despite 
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Fig. 6. This 12 CO (J=l-0) map dDame et alj|200lb shows the 
molecular clouds present around the coordinates of the burst 
(marked by the cross). Both the burst position and the distance 
to the dust-scattering layer suggest id entification of the laye r 
with the molecular cloud [KOY98] 66 (IKawamura et al . 1998). 
The black line marks the position of the galactic plane. 



GRB 070129 Galactic position 
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Fig. 7. A v map dSchlegel et alJI 1998b around the GRB 070129 
position (indicated by the cross). 



its smaller collecting area, Swift can also detect GRB halos, 
thanks to its capability to point in the GRB direction within a 
few tens of seconds, thus c atching the possible rin gs near the 
maximum of their intensity (IMereghetti et al.1 2006). 

Besides providing in some cases useful information on 
the GRB prompt X-ray emission, detection of these expand- 
ing dust rings offers a powerful tool for mapping the three- 
dimensional distribution of Galactic dust clouds. It should be 
noted that the derived distances, based purely on dynamical ar- 
guments, can be extremely accurate. 

This is exemplified by our results for GRB 061019. The 
line of sight of the burst passes through the Perseus spiral 
arm and near the galactic anti-center, where the usual kine- 
matic methods of estimating distances are not applicable. The 
molecular clouds and nebulae present in this region (see Figj6]l 



have been extensively studi ed (Lvnds 1962t iDame et al.ll2001 
lLeef l999; Ka wamura et alJll998l) . In particular, the position of 
GRB 061019 is consistent with the coordin ates of the cloud 
IKOY98] 66 of the 13 CO (J=l-0) survey by Kawamura et al. 
(119981) . placed at a distance of ~ 1.1 kpc from Earth. Thus our 



determination of the distance to the dust layer (D = 941^ pc) 
suggests an identification with the molecular cloud [KOY98] 
66, confirming and refining its distance estimate. 
Assuming dust properties, one can extract some information 
about the prompt X-ray emission of the GRB from the spec- 
trum and intensity of the halo, which is generally difficult 
to observe directly. This is exemplified by our study of the 
GRB 061019 halo: adopting the optical depth-energy rela- 
tion in Eq. [2] we found a prompt X-ray spectrum consistent 
with the extrapolation to the few keV rang e of the prompt 



hard X -ray emission observed by Swift-BAT (ISakamoto et al 
2006). Note that this apparently unsurprising r esult is not the 



rule, as shown by the case of G RB 031203 dVaughan et al 



2004; iTiengo & Mereghettil l2006). for which, independent of 



the adopted dust models, the derived soft X-ray emission is 
not compatible with an extrapolation of the hard X-ray spec- 
trum. 

When the prompt X-ray emission is known, the study of the 
rings in principle permits us to verify and eventually dis- 
tinguish between different proposed scattering optical depth- 
energy relations (like Eqsf2]and[3]), which depends via scatter- 
ing cross section on the dust properties like size and composi- 
tion. This is exemplified in our study of the rings around GRB 
070129, which could not however distinguish between the two 
proposed formula, due to the inadequate statistics of this halo. 
Future observations of prompt X-ray emissions followed by 
dust scattering rings as intense as those of GRB 061019 or GRB 
031203 might lead to more concrete possibilities of obtaining 
useful information on the properties of the interstellar dust in 
our Galaxy. 
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